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lamellar phases. 
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I. INTRODUCTION 


In order to reach a thorough understanding of colloidal suspensions, it is crucial to probe 
their dynamics at length scales comparable to the inter-particle distance. Thus, a consid¬ 
erable amount of experimental and theoretical work has been concerned with the dynamics 
of dense colloidal suspensions^. Space- and time-resolved experiments were mainly per¬ 
formed using dynamic light scattering (DLS) on micron-size spherical particles, both in the 
isotropio^ and the crystalline^ phases. 

DLS was also used to study dilute isotropic suspensions of elongated particles^ or solutions 
of disk-shaped particles close to the sol-gel transition!^, but using this technique for the 
study of concentrated and/or ordered solutions becomes challenging, due to the required 
particle sizes. This is nevertheless a regime where the system should exhibit rich dynamics, 
influenced by the (short- or long-range) order and by the hydrodynamic coupling, the latter 
being quantihed by the hydrodynamic function. 

Theoretical efforts aiming to calculate this function have focussed mainly on suspensions 
of spherical particles^; they have been largely validated by the experimental results (see 
below). For anisotropic particles, on the other hand, the calculations are much more in¬ 
volved and they were mostly restricted to solutions of slender rods^*^. We are not aware of 
any analytical results for the hydrodynamic function of plate-like particles. Experimental 
measurements at high scattering vector (corresponding to the interparticle distance) are also 
lacking, due to the intrinsic wavelength limitation of DLS. 

In recent years, the accessible distance range increased signihcantly via the use of X-ray 
Photon Correlation Spectroscopy (XPCS), which is fundamentally similar to DLS but uses 
X-ray radiation as a probe. Although this technique is only effective on systems with high 
scattering contrast and slow relaxation rates, it has already been used to study concentrated 
hard sphere solutions^i^— , aging suspensionsl^il^, particles dispersed in complex fluids^'^, 
and interface dynamics^il^. 

Some XPCS studies were performed on mesophases. For instance, the collective diffusion 
coefficient of a colloidal nematic phase was determined over a wide range of the wave vector 
q. The analysis showed that the dynamics of the system displays strong slowing down over 
length scales larger than the interparticle distancel^. The relaxation of capillary surface 
waves has also been measuredl^. Smectic phases -which have a certain degree of positional 


2 


order- represent ideal candidates for XPCS investigations, as the scattering intensity is 
concentrated in the vicinity of the Bragg peaks, leading to high count rates. That is why 
some of the hrst XPCS results were obtained on such systems, in particular under the form 
of free-standing22i^— dims. The bulk dynamics of a lamellar lyotropic phase of surfactant 
were also studied^^. The results, conhrmed by DLS measurements, could be interpreted in 
terms of the continuum elastic theory of smectics^iiS^. 

In this paper we present results obtained on lamellar lyotropic phases composed of large 
inorganic colloids; in this system, the slow relaxation rates and the high scattering contrast 
greatly extend the accessible g-range. We measure the collective diffusion coefficient D{qz), 
where i is the direction of the lamellar director, for a fairly wide range of scattering vectors 
along the lamellar director, (covering at least the hrst Bragg peak of the phase). We 
then obtain the hydrodynamic function H{qz). The large aspect ratio of the particles and 
their lamellar order enable us to describe the dynamics using a simplihed analytical model, 
which is in semi-quantitative agreement with the data. As a function of the concentration, 
the hydrodynamic coupling goes from very weak to extremely strong. 


II. EXPERIMENTAL 


We used concentrated suspensions of phosphatoantimonate (H3Sb3P20i4) particles, with 
a typical thickness of 1 nm and at least 300 nm wide^, synthesized as in reference l27|. The 
lateral size distribution of the particles was assessed using a scanning electron microscope 
equipped with a held emission gun (SEM-FEG). In a typical experiment, a colloidal sus¬ 
pension was diluted down to a volume fraction of 2.5 10“^ by adding distilled water. Then, 
a drop of suspension was spin-coated at ~ 1000 rpm onto a silicon wafer. Close inspection 
of the samples, by SEM-FEG at 10 kV and by optical microscopy, revealed the presence of 
plate-like particles that sometimes had a clear hexagonal shape. Their size distribution is 
extremely broad and ranges from a few hundred nanometers to more than hve microns. 

The batch solutions were diluted, either with pure water or with a 50:50 (wt%) glyc¬ 
erol/water mixture to reach the desired volume concentration 0. The room-temperature 
phase diagram was determined by visual inspection of the vials, in natural light and between 
crossed polarizers. The birefringent lamellar phase has a gel-like texture above 0 = 1.8% 
(does not how under its own weight, i.e. it has a zero-frequency elastic modulus) and is huid 
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below this value, down to 0 = 0.65% where it coexists with a very dilute isotropic phase, 


rase boundaries are compatible with the phase 


26l. for the case of no added salt (ionic strength 


with a concentration 0/ < 0.1%. These p 
diagram presented in Figure 3 of reference 
below 10“^ mol/1). 

The lamellar samples were prepared by aspiration in round glass capillaries, 1 mm wide 
and with a wall thickness of 10 /rm. By scanning the capillaries in the beam, it is easy to 
hnd aligned domains (usually, with the director normal to the capillary walls). 

The small-angle X-ray scattering (SAXS) and XPCS measurements were performed at 
the TROIKA beam line IDIOA of the ESRF with an X-ray energy of 8 keV (A = 1.55 A) 
selected by a single-bounce Si(lll) monochromator, in the uniform Filing mode of the storage 
ring. A (partially) coherent beam is obtained by inserting a 10 /rm pinhole aperture a few 
centimeters upstream of the sample. 

We used a 2D Maxipix detector consisting of 256 x 256 square pixels (55 pm in size) 
and with acquisition rates of up to 1 kHz. For XPCS, a few thousand images were recorded 
and the intensity autocorrelation functions were calculated by ensemble averaging^^ over 
pixels with the same absolute value of the scattering vector (see Figure [1]). The SAXS 
pattern is obtained as the average intensity over all the images in each run. Since we are 
only concerned with the Qz dependence of the measured quantities, in the following we will 
denote throughout by q. 

The relaxation of concentration fluctuations with a given wave vector q is reflected in the 
field correlation function at that value of the scattering vector: gi{q,t) = {E*{0)E{t)) / {E*{0)E{0)). 
As we will see below, our data is well described by a single exponential relaxation: 
gi{q,t) = exp[-D(g)t]. 

Experimentally, we measure the normalized intensity correlation function g 2 {q,t), related 
to the held correlation function by the Siegert relation: 


92 {q,t) = l + b{q) \gi{q,t)f = 1 + 6(g) exp[-2D(g)f] 


( 1 ) 


where the contrast factor 6 is a few percent. Fitting the experimental data to Equation 
m, as illustrated in Figure [H yields the relaxation rates D(g). We will further dehne a 
scale-dependent diffusion constant, D{q) = D(g)/g^. In the following, we will work with the 
(static) structure factor S{q) and with the diffusion constant D{q). For the sample with 
0 = 2.1 vol% (used as an example in Figures [H |2]and[3]) these parameters are displayed in 
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pixels 

FIG. 1. SAXS scattering pattern for the sample with cf) = 2.1 vol%; the lamellar director is roughly 
horizontal. The feature at the center of the image is the first Bragg sheet. The shadow of the 
beamstop is visible at the left edge. The averaging was done within the range delimited by the 
two dashed lines, on circular arcs corresponding to a given scattering vector q. One arc is shown 
as solid line, for q = 14.5 10“^ A ^. 

Figure |3l The diffusion constant is shown for all 0 values in Figure 01 


III. RESULTS AND DISCUSSION 


The structure factor and fluctuations of the smectic phase have been studied in detail over 


the last decades by many authors (see Ref. 


29 


for an in-depth review). The deformation free 


energy of the phase consists essentially of two terms, proportional to the compression modu¬ 
lus (along the smectic director) and to the bending stiffness of the layers. Both the amplitude 
and the relaxation of the deformation eigenmodes (corresponding to a unique phonon) can 
be described in a relatively straightforward (although algebraically cumbersome) fashion. 

However, if one is interested in the dynamic structure factor, i.e. in the collective behavior 
of the system at a given wave vector, the calculations become extremely involved, since all 
the phonons contribute in non-trivial ways. In order to make the analytical treatment as easy 
as possible, so that the underlying physics is not obscured by the mathematical formalism. 
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FIG. 2. Normalized intensity correlation function g 2 iq,t) for the sample with cj) = 2.1 vol% (see 
Figured]), shifted vertically in steps of 0.3. The different symbols correspond to different values of 
the scattering vector q. The solid lines are fits to Equation ([T]). 


some simplifying assumptions must be made. 


In this work, we are dealing with rather dilute phases composed of very stiff and very 
large platelets. We will therefore assume that the compression modulus is low and the 
bending stiffness is large, so that the most important fluctuations are those affecting the 
spacing between platelets along the director z and we treat the system as effectively one¬ 
dimensional, within the framework of a damped harmonic chain model. 
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FIG. 3. Structure factor S{q) and diffusion constant D{q) for the sample with cj) = 2.1 vol%. D^o 
(materialized by a dotted line) is the high-g value of the diffusion constant, corresponding to the 
range where the structure factor saturates to 1. 


A. Statics 


For this model, the static structure factor has an analytical solution, given by Refs, 
and l3ll : 

sinh ^ 


30 


/(g,0) = ^(g) = 


9 O' 
2 


cosh j — cos (qd) 


( 2 ) 


with a the typical fluctuation amplitude, dehned by = ksT /a as a function of the strength 
of the harmonic potential, a. 

Clearly, the harmonic model is not very accurate at high q, as one can see in Figure 
|5] for the hts to the low-concentration data. The interaction potential is probably stiffer, 
leading to pronounced second- and even third-order peaks, which are better described by 
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FIG. 4. Scale-dependent diffusion constant D{q) for different cj) values. The “dips” visible in the 
curves with cp = 0.8 and 1.1 vol% (and, to a lesser degree, for cp = 1.4 vol%) occur at the position 
of the peak in the static structure factor S{q). 

a hard-particle model^ (plotted as dashed lines). In the following we will nevertheless use 
the harmonic approximation, which describes the first peak rather well and which is more 
tractable as far as the dynamics are concerned (see Section fill Bj) . 

From the analysis of the structure factor data in Figure [5] using the model ([2]) we extract 
the repeat distance d and the fluctuation amplitude a. In Figure [6] we present both the 
dilution law d{(p) and the ratio cr/d as a function of the volume fraction (p. Three important 
conclusions can be drawn from this data: 

• The platelet thickness 6, given by the slope of the dilution law, is 7.4 ± 1.2A, much 
smaller than the 10 A value given in the literature^^. This discrepancy could be due 
to imperfect exfoliation of the layers, leading to a lower effective concentration in the 
lamellar phase. 

• There is no systematic difference in d and a between the systems formulated in pure 
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FIG. 5. Structure factor S{q) for different concentrations cj) of the lamellar phase (diamonds). 
Curves shifted vertically in steps of 1. The solid lines are fits to the harmonic chain model ([2]). For 
the two lowest concentrations, we also present as dashed lines the fit with the hard-rod model (see 
the text for more details). 


water (at room temperature) and those with 50 % glycerol (at 253 K). Since the phase 
boundaries are also the same, we conclude that the static structure of the phase is 
not affected by the presence of glycerol. At the same time, there is no discernible 
difference between the two synthesis batches. 


• The ratio a/d = 0.2 ±0.05 remains almost constant as (j) (and hence d) vary by almost 
an order of magnitude. 
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FIG. 6. Top: dilution law. The smectic repeat distance d is plotted against the volume fraction 
4> in log-log representation. The solid line represents the d ~ 0“^ variation expected for a one¬ 
dimensional system. The resulting platelet thickness, 7.4 ± 1.2 A, is clearly below the 10 A value 
given in the literature^^. Bottom: the ratio between the fluctuation amplitude a (see Equation [2|) 
and d. The data is shown for all the available samples, both in pure water at 293 K and in a 50:50 
glycerol/water mixture at 253 K. 


B. Dynamics 

To investigate the dynamics of the system we consider a damped harmonic chain, following 
the notations of Geisel^^. The chain consists of N particles at positions Xj, and we dehne the 
displacements Uj = Xj — x° with respect to the reference positions, which obey — Xj=d, 
where d is the average interparticle distance. For dehniteness, we assume periodic boundary 
conditions: un+i = Ui and only consider odd values for N. 

We start by considering the case of hydrodynamically uncoupled particles, i.e. the energy 
dissipation for each particle only depends on its velocity with respect to the surrounding 
fluid, and not on the position or velocity of the other particles. 
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The equations of motion are: 


dV , 

mxj -|- m'yxj + —— = r At) 
dxj 

(3a) 

1 ^ = ^ ^ {xk+i -Xk- df 

(3b) 

(Fj(t)Fi(O)) = 2m'ykBT6ijS{t) 

(3c) 


where m is the particle mass, mq is a damping coefficient (to be discussed in detail further 
on) and Fj{t) is the random force acting on particle j at time t. 

The quantity of interest is the time-dependent correlation function (depending also on 
the wave vector q): 

1 ^ 

/(g, t) = — J2{ (4) 

' l,k=l 

which reduces for t = 0 to the static structure factor (|2]). 

The treatment is quite standard: the displacements Uj are developed over the basis 
of normal modes (phonons). The amplitudes and time relaxation rates of the phonons 
obtained from (E]) are used to express I{q,t) by developing the right-hand side in Eq. (|lj). 
In the following, we restrict ourselves to the overdamped case, 7 S> ^/aJm = Uq. In 
this limit, the inertial term in Eq. fl3a| can be neglected and a unique relaxation rate is 
associated to each phonon^^. The resulting expression for /(g, t) is unwieldy, so we use the 
cumulant expansiort^^. The first cumulant, dehned as Ki{q) = 1 i-o represents 

an average relaxation rate for a given wave vector. Our experimental data is well described 
by a single exponential decay, I{q,t) = /(g, 0) in which case the first cumulant is 

just the relaxation rate: Ki{q) = fl{q), an assumption we will make throughout the analysis 
below. 

In the uncoupled case, one has simply: 


Kdl) 


with 




S{q) 


0 2 

q 


Dn = a 


2^0 


ksT 


(5) 


7 m7 

The relaxation rate obeys a diffusive law, with a diffusion constant D{q) = Do/S{q) 
which contains explicitly the well-known 1/S{q) dependence (“de Gennes narrowing”—). 
Dq is the value of D{q) at high wave vectors, where S{q) reaches 1, so that one should have 
D{q) ■ S(q)/D q = 1 at all q. In other words, only the thermodynamic forces (related to S{q)) 
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are taken into account, and the prefactor Dq is simply proportional to the mobility of an 
isolated particle, For a thin circular plate of radius R moving normal to its plane 

at low Reynolds number in a fluid with viscosity y, one has^: my = IQrjR. 



FIG. 7. Experimental data for the rescaled hydrodynamic function H{q)/H{oo) (symbols) and the 
static structure factor S{q) (solid lines) at various concentrations cj) indicated alongside the curves. 
The left and right panel are plotted against a common y-axis. For clarity, the curves are shifted 
vertically in steps of 3 units. 

However, one must also account for the hydrodynamic interactions, which lead to a more 
complex, and generally scale-dependent, form for the dissipation. Their effect is quantihed 
in terms of the hydrodynamic function, H{q), dehned as^: 

H(q) = (6) 

-^0 

In dilute solutions, H{q) ^ 1 at high q values (where S{q) saturates at 1) but is in general 
different from unity below this range. We emphasize that in the simplihed model discussed 
above, H{q) = 1 identically (from equations ([5]) and ([6]).) 
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Hydrodynamic interactions are essential in the dynamics of colloidal suspensions, even 
at moderate concentrations. In three-dimensional systems they have been studied quite 
extensively, both for spherical^ and for elongated^>^ particles. In recent years, the two- 
and one-dimensional cases have also been considered, in particular in the context of con- 
hned spherical particles^'—. In contrast, the phase under investigation consists of very 
anisotropic platelets and the one-dimensional character is given by the orientational order 
of the phase, rather than by conhnement; hence, the models cited above are probably not 
adapted. Since the distance between platelets along the smectic director (face to face) is 
much lower than their lateral extension, we only consider hydrodynamic coupling along this 
direction, neglecting the effect of the in-plane neighbors (edge to edge). 

The simplest model that accounts for these features is the well-known Stefan equation:^ 
describing the squeezing flow between two parallel plates at low Reynolds number (in the 
lubrication approximation). The viscous force acting on the plates is: 

37Tr]R^h 

with Tj the viscosity of the fluid, R the radius of the (circular) plates and h the gap width. 

In this case, the dissipation is no longer proportional to the velocity of an individual 
particle with respect to the surrounding fluid, but rather to the velocity difference between 
neighboring particles: 


T7V1SC 

= 


where 


-C + Xj+i)] 

^TirfR^ 


( 8 ) 


For this reason, the dissipation associated to the homogeneous (g = 0) mode vanishes: the 
particles can move together, at any velocity and maintaining the spacing d. The consequence 
is an inhnite relaxation rate. To remove this artihcial divergence, we also preserve an 
individual friction term (formally identical to the second term on the Ihs in Equation fl3a)) .) 
In this case, however, the quantitative value for my is different from the free case, since the 
interaction of each particle with the fluid is “screened” by its neighbors. A simple yet fairly 
realistic model is that of a cylindrical stack of platelets: C corresponds then to the relative 
motion of particles within the stack, while my is associated with a solid-like translation of 
the cylinder along its axis, the dissipation taking place in the surrounding medium (with 
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an effective viscosity rjes which is in general different from that of the solvent). Using well- 
known formnlas for the longitndinal mobility of a cylinder— and considering an aspect ratio 
(length L to diameter 2R) of the order of 10, one has: D\\ = amonnting to a friction 

coefficient per particle my = 

With the ingredients above, and after some algebra (detailed in the Appendix) one obtains 
for H{q) an explicit formula: 


Af-l, 


H{q)= 




(l - cos(gdp)exp 




P X 


xz^ , with 2 ; = 1 + 1/13 — a/1 -|- 2/3//9 


(9) 


where N is the number of particles and (3 = 2(^/[rri'y) is the ratio of the collective dissipation 
described by ([8]) to the individual particle dissipation; it provides a quantitative measure of 
the hydrodynamic coupling. The primed sum symbol indicates that the p = 0 term should 
be halved. In the high-g limit, H{q) — H(oo) = {2f3 + 1 )“^A so that the corresponding 
diffusion constant is: 


^ _ -Do _ ksT 

°° \/2/3 -I- 1 m'y4:(/{m'y) + 1 


( 10 ) 


As shown in the Appendix, for high values of the coupling constant (3 the rescaled hydro- 
dynamic function H{q)/H{oo) approaches the static structure factor S'(g); in other words, 
D(g) becomes constant since the hydrodynamic effects are much stronger than the ther¬ 
modynamic ones. This tendency is indeed visible in Figure IH where the “dips” present at 
the peak positions for low concentrations flatten out^. We would then expect that, as the 
concentration increases, the hydrodynamic function evolves from a constant value to the 
limiting shape S{q). This is indeed observed in Figure [71 where there is a clear difference 
between the data at 0 = 0.8 % and those at higher concentration. Indeed, at cf) = 1.1 % 
H{q) is already similar to S{q) and the similarity becomes clearer above this 0 value. More 
points within 0.8 % < 0 < 1.1 % would be needed to resolve the onset of coupling. 

The similarity between S{q) and H{q) is very useful, since the static structure factor can 
be measured much more precisely and to higher q values than the hydrodynamic function. 
We therefore rescale H{q), bringing it as close as possible to S{q) (the data in Figure [7] 
has already been rescaled). This operation also sets the value of Dqo without the need of 
measuring D{q) to very high q. 

To summarize, the high-g value of the diffusion coefficient in the two regimes is given by: 
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FIG. 8. High-g value of the diffusion constant, D^o, as a function of the volume fraction cj) of 
platelets in a 50:50 glycerol/water mixture at 253 K. The phase diagram is also indicated: the 
(fluid or gel-like) lamellar phase at high concentration and an isotropic (I) phase at very low 
concentration, with a rather wide coexistence range. 
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where we considered that the dissipation by squeezing flow between the particles involved 
the solvent viscosity ?7 = 50 mPa s (for a 50:50 wt % mixture of glycerol in water at 253 K— ) 
while for the dissipation of the individual particles one needs to use the effective viscosity 
of the medium, taken as independent of the concentration in the investigated range. 
The smectic repeat distance d is given in Figure El so the only adjustable parameters are 
r^eff and the platelet radius R. Good agreement with the experimental data is obtained with 
the values: rjes = 3.27] and R = 3.5/im (see Figure [H]). Indeed, at the lowest concentration 
0 = 0.8 % one has D^o = Dq given by Eq. flllap (dashed line in Figure |8]), while at higher 
concentration the data is well described by the dependence flllbp . plotted as solid line. The 
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coupling occurs for 0 between 0.8 and 1.1 %, in agreement with the interpretation proposed 
above for the hydrodynamic function. It is noteworthy that no signihcant change in the 
microscopic dynamics can be detected at the transition between the fluid and the gel-like 
lamellar phases. 

The rather large value of the platelet radius obtained above is not surprising, since the 
largest particles dominate both the scattering signal (with a contribution proportional to 
the square of the particle volume) and the correlation function, contributing the slowest 
relaxation. Furthermore, the coupling in Equation flllbp goes as the fourth power of the 
radius. The size distribution being very wide (see Section HI]), all experimental results should 
be severely skewed towards the large particles. 


IV. CONCLUSION 

We measured the static and dynamic properties of a lamellar phase composed of rigid 
platelets and quantihed the hydrodynamic coupling between nearest neighbors (along the 
director of the phase). The coupling is almost absent at low concentrations, where the 
dissipation occurs at the level of the isolated particle, but it quickly becomes dominant at 
higher concentrations. The hydrodynamic function of the phase H{qz) is relatively well 
described by an analytical one-dimensional model. 

This result is noteworthy insofar the hydrodynamic function -although indispensable for 
understanding the relaxation at the particle size- is generally difficult to calculate, even for 
spherical colloids. To our knowledge, no explicit models have been proposed for anisotropic 
particles. Somewhat surprisingly, the order of the particles (which should further complicate 
the study) allows us in the present case to use a very simplified approach. 

The high-g value of the diffusion coefficient Doo decreases rapidly with the concentration 
in the fluid lamellar phase, but it remains practically constant across the sol/gel transition 
and well into the gel-like regime. No sharp variation is observed at the transition, in either 
the dynamic or the static parameters (in particular, we do not detect any spatial inho¬ 
mogeneities), showing that this transition involves longer length- and/or time-scales than 
probed in our experiment. 
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Appendix: Coupled hydrodynamics 

In this Appendix we detail the calcnlations of the hydrodynamic fnnction for the harmonic 
model (jS]), with or withont the conpling term (|8]). 

The bonndary conditions specihed in section IIII Bl impose the eigenvector basis: 

271 n A — 1 , , ,, /. X 

Qn = n = 0, ±1,... ± ——, (A odd) (A.l) 

snch that —7i/d < qn < vr/d (restriction to the hrst Brillonin zone). The normal modes 
(phonons) are given by fn{xi) = exp{iqndl) = exp{2i7iln/N) and the individnal displace¬ 
ments are expressed as: uft) = Ylin^riif)fn{.xi)- 
The amplitnde coefficients obey: 




knT 


4Aa sin^(q'„d/2) 

Normal mode expansion of Eqnation fl3al) with the additional dissipative term 
the relaxation rate of the phonons: 


(A.2) 

yields 
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{A*^{Q)An{t)) = (5m,n (|^n|^)exp(-r„t), 
2a 1 - cos{qnd) 

with i n = 


(A.3) 


rri'-y 1 + /3 [1 — cos(g„(i)] 

where we remind that (3 = 2C,/{rri'y) is the ratio of the collective dissipation to the individual 
particle dissipation. 

The time-dependent correlation function is then obtained by simple substitution in (jTj): 

N-l, 


I{q, t)= 2 ^ ^ cos{qdp) exp 


X 


p=0 


-V 


N 


q^a^ 


■X 


1 - cos{qndp) exp{-rnt) 


N ^ 1 — cos{qnd) 

and the hrst cumulant is obtained according to ([3]), yielding for H{q): 


(A.4) 


^ N — p f q'^a'^ 

H(q)= E 2 — cos{qdp) exp (- —p ] x 


p=0 

xC(p,/3) 


N 


(A.5) 


where the coefficients C{p,f3) = i+/ 3 [r-cosfg d)] reduced (e.g. by going to the 

continuum limit and performing a complex integral over the unit circle) to: C{p, (3) = 
finally yielding equation (|9]) above. For the uncoupled case /9 = 0 this result is greatly 
simplified, since the sum over the modes in Equation ([9]) becomes: ^ Yhn cos(g„(ip) = 
and thus H{q, /3 = 0) = 1, leading to the second equality in Equation (j5]). 

Q / \ 

The strongly coupled form lim H(q) = - is obtained by noting that in this limit 

/3^oo v 2/3 -I-1 


-{^-\nz)p 


is dominated 


z approaches 1 from below, so that the damping factor exp 
by the first term in the exponent. One can therefore drop the factor in Eq. (1A.5p and 
retrieve an obvious expansion of the static structure factor S{q). Of course, for this to 
hold the scattering vector must also be above a certain value, otherwise the first term in the 
exponent might be too small. We checked that for reasonable numerical values the condition 
is fulfilled for all accessible values of q. 
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